Expanding the repertoire reveals recurrent, cryptic, and hematopoietic HLA class I minor histocompatibility antigens.

ABSTRACT: Allogeneic stem cell transplantation (alloSCT) is a curative treatment for hematological malignancies. After HLA-matched alloSCT, antitumor immunity is caused by donor T cells recognizing polymorphic peptides, designated minor histocompatibility antigens (MiHAs), that are presented by HLA on malignant patient cells. However, T cells often target MiHAs on healthy nonhematopoietic tissues of patients, thereby inducing side effects known as graft-versus-host disease. Here, we aimed to identify the dominant repertoire of HLA-I-restricted MiHAs to enable strategies to predict, monitor or modulate immune responses after alloSCT. To systematically identify novel MiHAs by genome-wide association screening, T-cell clones were isolated from 39 transplanted patients and tested for reactivity against 191 Epstein-Barr virus transformed B cell lines of the 1000 Genomes Project. By discovering 81 new MiHAs, we more than doubled the antigen repertoire to 159 MiHAs and demonstrated that, despite many genetic differences between patients and donors, often the same MiHAs are targeted in multiple patients. Furthermore, we showed that one quarter of the antigens are cryptic, that is translated from unconventional open reading frames, for example long noncoding RNAs, showing that these antigen types are relevant targets in natural immune responses. Finally, using single cell RNA-seq data, we analyzed tissue expression of MiHA-encoding genes to explore their potential role in clinical outcome, and characterized 11 new hematopoietic-restricted MiHAs as potential targets for immunotherapy. In conclusion, we expanded the repertoire of HLA-I-restricted MiHAs and identified recurrent, cryptic and hematopoietic-restricted antigens, which are fundamental to predict, follow or manipulate immune responses to improve clinical outcome after alloSCT.

PANDORA v2.0: Benchmarking peptide-MHC II models and software improvements.

T-cell specificity to differentiate between self and non-self relies on T-cell receptor (TCR) recognition of peptides presented by the Major Histocompatibility Complex (MHC). Investigations into the three-dimensional (3D) structures of peptide:MHC (pMHC) complexes have provided valuable insights of MHC functions. Given the limited availability of experimental pMHC structures and considerable diversity of peptides and MHC alleles, it calls for the development of efficient and reliable computational approaches for modeling pMHC structures. Here we present an update of PANDORA and the systematic evaluation of its performance in modelling 3D structures of pMHC class II complexes (pMHC-II), which play a key role in the cancer immune response. PANDORA is a modelling software that can build low-energy models in a few minutes by restraining peptide residues inside the MHC-II binding groove. We benchmarked PANDORA on 136 experimentally determined pMHC-II structures covering 44 unique αß chain pairs. Our pipeline achieves a median backbone Ligand-Root Mean Squared Deviation (L-RMSD) of 0.42 Å on the binding core and 0.88 Å on the whole peptide for the benchmark dataset. We incorporated software improvements to make PANDORA a pan-allele framework and improved the user interface and software quality. Its computational efficiency allows enriching the wealth of pMHC binding affinity and mass spectrometry data with 3D models. These models can be used as a starting point for molecular dynamics simulations or structure-boosted deep learning algorithms to identify MHC-binding peptides. PANDORA is available as a Python package through Conda or as a source installation at https://github.com/X-lab-3D/PANDORA.

PANDORA: A Fast, Anchor-Restrained Modelling Protocol for Peptide: MHC Complexes.

Deeper understanding of T-cell-mediated adaptive immune responses is important for the design of cancer immunotherapies and antiviral vaccines against pandemic outbreaks. T-cells are activated when they recognize foreign peptides that are presented on the cell surface by Major Histocompatibility Complexes (MHC), forming peptide:MHC (pMHC) complexes. 3D structures of pMHC complexes provide fundamental insight into T-cell recognition mechanism and aids immunotherapy design. High MHC and peptide diversities necessitate efficient computational modelling to enable whole proteome structural analysis. We developed PANDORA, a generic modelling pipeline for pMHC class I and II (pMHC-I and pMHC-II), and present its performance on pMHC-I here. Given a query, PANDORA searches for structural templates in its extensive database and then applies anchor restraints to the modelling process. This restrained energy minimization ensures one of the fastest pMHC modelling pipelines so far. On a set of 835 pMHC-I complexes over 78 MHC types, PANDORA generated models with a median RMSD of 0.70 Å and achieved a 93% success rate in top 10 models. PANDORA performs competitively with three pMHC-I modelling state-of-the-art approaches and outperforms AlphaFold2 in terms of accuracy while being superior to it in speed. PANDORA is a modularized and user-configurable python package with easy installation. We envision PANDORA to fuel deep learning algorithms with large-scale high-quality 3D models to tackle long-standing immunology challenges.

Genome-wide identification of genes regulating DNA methylation using genetic anchors for causal inference.

BACKGROUND: DNA methylation is a key epigenetic modification in human development and disease, yet there is limited understanding of its highly coordinated regulation. Here, we identify 818 genes that affect DNA methylation patterns in blood using large-scale population genomics data. RESULTS: By employing genetic instruments as causal anchors, we establish directed associations between gene expression and distant DNA methylation levels, while ensuring specificity of the associations by correcting for linkage disequilibrium and pleiotropy among neighboring genes. The identified genes are enriched for transcription factors, of which many consistently increased or decreased DNA methylation levels at multiple CpG sites. In addition, we show that a substantial number of transcription factors affected DNA methylation at their experimentally determined binding sites. We also observe genes encoding proteins with heterogenous functions that have widespread effects on DNA methylation, e.g., NFKBIE, CDCA7(L), and NLRC5, and for several examples, we suggest plausible mechanisms underlying their effect on DNA methylation. CONCLUSION: We report hundreds of genes that affect DNA methylation and provide key insights in the principles underlying epigenetic regulation.

Tumor Necrosis Factor Receptor SF10A (TNFRSF10A) SNPs Correlate With Corticosteroid Response in Duchenne Muscular Dystrophy.

BACKGROUND: Duchenne muscular dystrophy (DMD) is a rare and severe X-linked muscular dystrophy in which the standard of care with variable outcome, also due to different drug response, is chronic off-label treatment with corticosteroids (CS). In order to search for SNP biomarkers for corticosteroid responsiveness, we genotyped variants across 205 DMD-related genes in patients with differential response to steroid treatment. METHODS AND FINDINGS: We enrolled a total of 228 DMD patients with identified dystrophin mutations, 78 of these patients have been under corticosteroid treatment for at least 5 years. DMD patients were defined as high responders (HR) if they had maintained the ability to walk after 15 years of age and low responders (LR) for those who had lost ambulation before the age of 10 despite corticosteroid therapy. Based on interactome mapping, we prioritized 205 genes and sequenced them in 21 DMD patients (discovery cohort or DiC = 21). We identified 43 SNPs that discriminate between HR and LR. Discriminant Analysis of Principal Components (DAPC) prioritized 2 response-associated SNPs in the TNFRSF10A gene. Validation of this genotype was done in two additional larger cohorts composed of 46 DMD patients on corticosteroid therapy (validation cohorts or VaC1), and 150 non ambulant DMD patients and never treated with corticosteroids (VaC2). SNP analysis in all validation cohorts (N = 207) showed that the CT haplotype is significantly associated with HR DMDs confirming the discovery results. CONCLUSION: We have shown that TNFRSF10A CT haplotype correlates with corticosteroid response in DMD patients and propose it as an exploratory CS response biomarker.

Optimized Whole Genome Association Scanning for Discovery of HLA Class I-Restricted Minor Histocompatibility Antigens.

Patients undergoing allogeneic stem cell transplantation as treatment for hematological diseases face the risk of Graft-versus-Host Disease as well as relapse. Graft-versus-Host Disease and the favorable Graft-versus-Leukemia effect are mediated by donor T cells recognizing polymorphic peptides, which are presented on the cell surface by HLA molecules and result from single nucleotide polymorphism alleles that are disparate between patient and donor. Identification of polymorphic HLA-binding peptides, designated minor histocompatibility antigens, has been a laborious procedure, and the number and scope for broad clinical use of these antigens therefore remain limited. Here, we present an optimized whole genome association approach for discovery of HLA class I minor histocompatibility antigens. T cell clones isolated from patients who responded to donor lymphocyte infusions after HLA-matched allogeneic stem cell transplantation were tested against a panel of 191 EBV-transformed B cells, which have been sequenced by the 1000 Genomes Project and selected for expression of seven common HLA class I alleles (HLA-A∗01:01, A∗02:01, A∗03:01, B∗07:02, B∗08:01, C∗07:01, and C∗07:02). By including all polymorphisms with minor allele frequencies above 0.01, we demonstrated that the new approach allows direct discovery of minor histocompatibility antigens as exemplified by seven new antigens in eight different HLA class I alleles including one antigen in HLA-A∗24:02 and HLA-A∗23:01, for which the method has not been originally designed. Our new whole genome association strategy is expected to rapidly augment the repertoire of HLA class I-restricted minor histocompatibility antigens that will become available for donor selection and clinical use to predict, follow or manipulate Graft-versus-Leukemia effect and Graft-versus-Host Disease after allogeneic stem cell transplantation.

Prioritization of novel ADPKD drug candidates from disease-stage specific gene expression profiles.

BACKGROUND: Autosomal Dominant Polycystic Kidney Disease (ADPKD) is one of the most common causes of end-stage renal failure, caused by mutations in PKD1 or PKD2 genes. Tolvaptan, the only drug approved for ADPKD treatment, results in serious side-effects, warranting the need for novel drugs. METHODS: In this study, we applied RNA-sequencing of Pkd1cko mice at different disease stages, and with/without drug treatment to identify genes involved in ADPKD progression that were further used to identify novel drug candidates for ADPKD. We followed an integrative computational approach using a combination of gene expression profiling, bioinformatics and cheminformatics data. FINDINGS: We identified 1162 genes that had a normalized expression after treating the mice with drugs proven effective in preclinical models. Intersecting these genes with target affinity profiles for clinically-approved drugs in ChEMBL, resulted in the identification of 116 drugs targeting 29 proteins, of which several are previously linked to Polycystic Kidney Disease such as Rosiglitazone. Further testing the efficacy of six candidate drugs for inhibition of cyst swelling using a human 3D-cyst assay, revealed that three of the six had cyst-growth reducing effects with limited toxicity. INTERPRETATION: Our data further establishes drug repurposing as a robust drug discovery method, with three promising drug candidates identified for ADPKD treatment (Meclofenamic Acid, Gamolenic Acid and Birinapant). Our strategy that combines multiple-omics data, can be extended for ADPKD and other diseases in the future. FUNDING: European Union's Seventh Framework Program, Dutch Technology Foundation Stichting Technische Wetenschappen and the Dutch Kidney Foundation.

Characterisation of transcription factor profiles in polycystic kidney disease (PKD): identification and validation of STAT3 and RUNX1 in the injury/repair response and PKD progression.

Autosomal dominant polycystic kidney disease (ADPKD) is the most common genetic renal disease, caused in the majority of the cases by a mutation in either the PKD1 or the PKD2 gene. ADPKD is characterised by a progressive increase in the number and size of cysts, together with fibrosis and distortion of the renal architecture, over the years. This is accompanied by alterations in a complex network of signalling pathways. However, the underlying molecular mechanisms are not well characterised. Previously, we defined the PKD Signature, a set of genes typically dysregulated in PKD across different disease models from a meta-analysis of expression profiles. Given the importance of transcription factors (TFs) in modulating disease, we focused in this paper on characterising TFs from the PKD Signature. Our results revealed that out of the 1515 genes in the PKD Signature, 92 were TFs with altered expression in PKD, and 32 of those were also implicated in tissue injury/repair mechanisms. Validating the dysregulation of these TFs by qPCR in independent PKD and injury models largely confirmed these findings. STAT3 and RUNX1 displayed the strongest activation in cystic kidneys, as demonstrated by chromatin immunoprecipitation (ChIP) followed by qPCR. Using immunohistochemistry, we showed a dramatic increase of expression after renal injury in mice and cystic renal tissue of mice and humans. Our results suggest a role for STAT3 and RUNX1 and their downstream targets in the aetiology of ADPKD and indicate that the meta-analysis approach is a viable strategy for new target discovery in PKD. KEY MESSAGES: We identified a list of transcription factors (TFs) commonly dysregulated in ADPKD. Out of the 92 TFs identified in the PKD Signature, 35% are also involved in injury/repair processes. STAT3 and RUNX1 are the most significantly dysregulated TFs after injury and during PKD progression. STAT3 and RUNX1 activity is increased in cystic compared to non-cystic mouse kidneys. Increased expression of STAT3 and RUNX1 is observed in the nuclei of renal epithelial cells, also in human ADPKD samples.

Evaluation of commonly used analysis strategies for epigenome- and transcriptome-wide association studies through replication of large-scale population studies.

BACKGROUND: A large number of analysis strategies are available for DNA methylation (DNAm) array and RNA-seq datasets, but it is unclear which strategies are best to use. We compare commonly used strategies and report how they influence results in large cohort studies. RESULTS: We tested the associations of DNAm and RNA expression with age, BMI, and smoking in four different cohorts (n = ~ 2900). By comparing strategies against the base model on the number and percentage of replicated CpGs for DNAm analyses or genes for RNA-seq analyses in a leave-one-out cohort replication approach, we find the choice of the normalization method and statistical test does not strongly influence the results for DNAm array data. However, adjusting for cell counts or hidden confounders substantially decreases the number of replicated CpGs for age and increases the number of replicated CpGs for BMI and smoking. For RNA-seq data, the choice of the normalization method, gene expression inclusion threshold, and statistical test does not strongly influence the results. Including five principal components or excluding correction of technical covariates or cell counts decreases the number of replicated genes. CONCLUSIONS: Results were not influenced by the normalization method or statistical test. However, the correction method for cell counts, technical covariates, principal components, and/or hidden confounders does influence the results.

Search for Early Pancreatic Cancer Blood Biomarkers in Five European Prospective Population Biobanks Using Metabolomics.

Most patients with pancreatic cancer present with advanced disease and die within the first year after diagnosis. Predictive biomarkers that signal the presence of pancreatic cancer in an early stage are desperately needed. We aimed to identify new and validate previously found plasma metabolomic biomarkers associated with early stages of pancreatic cancer. Prediagnostic blood samples from individuals who were to receive a diagnosis of pancreatic cancer between 1 month and 17 years after sampling (N = 356) and age- and sex-matched controls (N = 887) were collected from five large population cohorts (HUNT2, HUNT3, FINRISK, Estonian Biobank, Rotterdam Study). We applied proton nuclear magnetic resonance-based metabolomics on the Nightingale platform. Logistic regression identified two interesting hits: glutamine (P = 0.011) and histidine (P = 0.012), with Westfall-Young family-wise error rate adjusted P values of 0.43 for both. Stratification in quintiles showed a 1.5-fold elevated risk for the lowest 20% of glutamine and a 2.2-fold increased risk for the lowest 20% of histidine. Stratification by time to diagnosis suggested glutamine to be involved in an earlier process (2 to 5 years before diagnosis), and histidine in a process closer to the actual onset (<2 years). Our data did not support the branched-chain amino acids identified earlier in several US cohorts as potential biomarkers for pancreatic cancer. Thus, although we identified glutamine and histidine as potential biomarkers of biological interest, our results imply that a study at this scale does not yield metabolomic biomarkers with sufficient predictive value to be clinically useful per se as prognostic biomarkers.

Skewed X-inactivation is common in the general female population.

X-inactivation is a well-established dosage compensation mechanism ensuring that X-chromosomal genes are expressed at comparable levels in males and females. Skewed X-inactivation is often explained by negative selection of one of the alleles. We demonstrate that imbalanced expression of the paternal and maternal X-chromosomes is common in the general population and that the random nature of the X-inactivation mechanism can be sufficient to explain the imbalance. To this end, we analyzed blood-derived RNA and whole-genome sequencing data from 79 female children and their parents from the Genome of the Netherlands project. We calculated the median ratio of the paternal over total counts at all X-chromosomal heterozygous single-nucleotide variants with coverage ≥10. We identified two individuals where the same X-chromosome was inactivated in all cells. Imbalanced expression of the two X-chromosomes (ratios ≤0.35 or ≥0.65) was observed in nearly 50% of the population. The empirically observed skewing is explained by a theoretical model where X-inactivation takes place in an embryonic stage in which eight cells give rise to the hematopoietic compartment. Genes escaping X-inactivation are expressed from both alleles and therefore demonstrate less skewing than inactivated genes. Using this characteristic, we identified three novel escapee genes (SSR4, REPS2, and SEPT6), but did not find support for many previously reported escapee genes in blood. Our collective data suggest that skewed X-inactivation is common in the general population. This may contribute to manifestation of symptoms in carriers of recessive X-linked disorders. We recommend that X-inactivation results should not be used lightly in the interpretation of X-linked variants.

Comparative transcriptomics of shear stress treated Pkd1-/- cells and pre-cystic kidneys reveals pathways involved in early polycystic kidney disease.

Mutations in the PKD1 or PKD2 genes are the cause of autosomal dominant polycystic kidney disease (ADPKD). The encoded proteins localize within the cell membrane and primary cilia and are proposed to be involved in mechanotransduction. Therefore, we evaluate shear stress dependent signaling in renal epithelial cells and the relevance for ADPKD. Using RNA sequencing and pathway analysis, we compared gene expression of in vitro shear stress treated Pkd1-/- renal epithelial cells and in vivo pre-cystic Pkd1del models. We show that shear stress alters the same signaling pathways in Pkd1-/- renal epithelial cells and Pkd1wt controls. However, expression of a number of genes was slightly more induced by shear stress in Pkd1-/- cells, suggesting that Pkd1 has the function to restrain shear regulated signaling instead of being a mechano-sensing activator. We also compared altered gene expression in Pkd1-/- cells during shear with in vivo transcriptome data of kidneys from Pkd1del mice at three early pre-cystic time-points. This revealed overlap of a limited number of differentially expressed genes. However, the overlap between cells and mice is much higher when looking at pathways and molecular processes, largely due to altered expression of paralogous genes. Several of the altered pathways in the in vitro and in vivo Pkd1del models are known to be implicated in ADPKD pathways, including PI3K-AKT, MAPK, Hippo, calcium, Wnt, and TGF-ß signaling. We hypothesize that increased activation of selected genes in renal epithelial cells early upon Pkd1 gene disruption may disturb the balance in signaling and may contribute to cyst formation.

Genome-wide identification of directed gene networks using large-scale population genomics data.

Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene-gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identified 49 genes as drivers of downstream transcriptional changes (Wald P < 7 × 10-10), among which transcription factors were overrepresented (Fisher's P = 3.3 × 10-7). Our analysis suggests new gene functions and targets, including for SENP7 (zinc-finger genes involved in retroviral repression) and BCL2A1 (target genes possibly involved in auditory dysfunction). Our work highlights the utility of population genomics data in deriving directed gene expression networks. A resource of trans-effects for all 6600 genes with a genetic instrument can be explored individually using a web-based browser.

Strap associates with Csde1 and affects expression of select Csde1-bound transcripts.

Erythropoiesis is regulated at many levels, including control of mRNA translation. Changing environmental conditions, such as hypoxia or the availability of nutrients and growth factors, require a rapid response enacted by the enhanced or repressed translation of existing transcripts. Cold shock domain protein e1 (Csde1/Unr) is an RNA-binding protein required for erythropoiesis and strongly upregulated in erythroblasts relative to other hematopoietic progenitors. The aim of this study is to identify the Csde1-containing protein complexes and investigate their role in post-transcriptional expression control of Csde1-bound transcripts. We show that Serine/Threonine kinase receptor-associated protein (Strap/Unrip), was the protein most strongly associated with Csde1 in erythroblasts. Strap is a WD40 protein involved in signaling and RNA splicing, but its role when associated with Csde1 is unknown. Reduced expression of Strap did not alter the pool of transcripts bound by Csde1. Instead, it altered the mRNA and/or protein expression of several Csde1-bound transcripts that encode for proteins essential for translational regulation during hypoxia, such as Hmbs, eIF4g3 and Pabpc4. Also affected by Strap knockdown were Vim, a Gata-1 target crucial for erythrocyte enucleation, and Elavl1, which stabilizes Gata-1 mRNA. The major cellular processes affected by both Csde1 and Strap were ribosome function and cell cycle control.

Brain Transcriptomic Analysis of Hereditary Cerebral Hemorrhage With Amyloidosis-Dutch Type.

Hereditary cerebral hemorrhage with amyloidosis-Dutch type (HCHWA-D) is an early onset hereditary form of cerebral amyloid angiopathy (CAA) caused by a point mutation resulting in an amino acid change (NP_000475.1:p.Glu693Gln) in the amyloid precursor protein (APP). Post-mortem frontal and occipital cortical brain tissue from nine patients and nine age-related controls was used for RNA sequencing to identify biological pathways affected in HCHWA-D. Although previous studies indicated that pathology is more severe in the occipital lobe in HCHWA-D compared to the frontal lobe, the current study showed similar changes in gene expression in frontal and occipital cortex and the two brain regions were pooled for further analysis. Significantly altered pathways were analyzed using gene set enrichment analysis (GSEA) on 2036 significantly differentially expressed genes. Main pathways over-represented by down-regulated genes were related to cellular aerobic respiration (including ATP synthesis and carbon metabolism) indicating a mitochondrial dysfunction. Principal up-regulated pathways were extracellular matrix (ECM)-receptor interaction and ECM proteoglycans in relation with an increase in the transforming growth factor beta (TGFß) signaling pathway. Comparison with the publicly available dataset from pre-symptomatic APP-E693Q transgenic mice identified overlap for the ECM-receptor interaction pathway, indicating that ECM modification is an early disease specific pathomechanism.

Immune stimuli shape the small non-coding transcriptome of extracellular vesicles released by dendritic cells.

The release and uptake of nano-sized extracellular vesicles (EV) is a highly conserved means of intercellular communication. The molecular composition of EV, and thereby their signaling function to target cells, is regulated by cellular activation and differentiation stimuli. EV are regarded as snapshots of cells and are, therefore, in the limelight as biomarkers for disease. Although research on EV-associated RNA has predominantly focused on microRNAs, the transcriptome of EV consists of multiple classes of small non-coding RNAs with potential gene-regulatory functions. It is not known whether environmental cues imposed on cells induce specific changes in a broad range of EV-associated RNA classes. Here, we investigated whether immune-activating or -suppressing stimuli imposed on primary dendritic cells affected the release of various small non-coding RNAs via EV. The small RNA transcriptomes of highly pure EV populations free from ribonucleoprotein particles were analyzed by RNA sequencing and RT-qPCR. Immune stimulus-specific changes were found in the miRNA, snoRNA, and Y-RNA content of EV from dendritic cells, whereas tRNA and snRNA levels were much less affected. Only part of the changes in EV-RNA content reflected changes in cellular RNA, which urges caution in interpreting EV as snapshots of cells. By comprehensive analysis of RNA obtained from highly purified EV, we demonstrate that multiple RNA classes contribute to genetic messages conveyed via EV. The identification of multiple RNA classes that display cell stimulation-dependent association with EV is the prelude to unraveling the function and biomarker potential of these EV-RNAs.

Full-length mRNA sequencing uncovers a widespread coupling between transcription initiation and mRNA processing.

BACKGROUND: The multifaceted control of gene expression requires tight coordination of regulatory mechanisms at transcriptional and post-transcriptional level. Here, we studied the interdependence of transcription initiation, splicing and polyadenylation events on single mRNA molecules by full-length mRNA sequencing. RESULTS: In MCF-7 breast cancer cells, we find 2700 genes with interdependent alternative transcription initiation, splicing and polyadenylation events, both in proximal and distant parts of mRNA molecules, including examples of coupling between transcription start sites and polyadenylation sites. The analysis of three human primary tissues (brain, heart and liver) reveals similar patterns of interdependency between transcription initiation and mRNA processing events. We predict thousands of novel open reading frames from full-length mRNA sequences and obtained evidence for their translation by shotgun proteomics. The mapping database rescues 358 previously unassigned peptides and improves the assignment of others. By recognizing sample-specific amino-acid changes and novel splicing patterns, full-length mRNA sequencing improves proteogenomics analysis of MCF-7 cells. CONCLUSIONS: Our findings demonstrate that our understanding of transcriptome complexity is far from complete and provides a basis to reveal largely unresolved mechanisms that coordinate transcription initiation and mRNA processing.

Csde1 binds transcripts involved in protein homeostasis and controls their expression in an erythroid cell line.

Expression of the RNA-binding protein Csde1 (Cold shock domain protein e1) is strongly upregulated during erythropoiesis compared to other hematopoietic lineages. Csde1 expression is impaired in the severe congenital anemia Diamond Blackfan Anemia (DBA), and reduced expression of Csde1 in healthy erythroblasts impaired their proliferation and differentiation. To investigate the cellular pathways controlled by Csde1 in erythropoiesis, we identified the transcripts that physically associate with Csde1 in erythroid cells. These mainly encoded proteins involved in ribogenesis, mRNA translation and protein degradation, but also proteins associated with the mitochondrial respiratory chain and mitosis. Crispr/Cas9-mediated deletion of the first cold shock domain of Csde1 affected RNA expression and/or protein expression of Csde1-bound transcripts. For instance, protein expression of Pabpc1 was enhanced while Pabpc1 mRNA expression was reduced indicating more efficient translation of Pabpc1 followed by negative feedback on mRNA stability. Overall, the effect of reduced Csde1 function on mRNA stability and translation of Csde1-bound transcripts was modest. Clones with complete loss of Csde1, however, could not be generated. We suggest that Csde1 is involved in feed-back control in protein homeostasis and that it dampens stochastic changes in mRNA expression.

Systematic genomic and translational efficiency studies of uveal melanoma.

To further our understanding of the somatic genetic basis of uveal melanoma, we sequenced the protein-coding regions of 52 primary tumors and 3 liver metastases together with paired normal DNA. Known recurrent mutations were identified in GNAQ, GNA11, BAP1, EIF1AX, and SF3B1. The role of mutated EIF1AX was tested using loss of function approaches including viability and translational efficiency assays. Knockdown of both wild type and mutant EIF1AX was lethal to uveal melanoma cells. We probed the function of N-terminal tail EIF1AX mutations by performing RNA sequencing of polysome-associated transcripts in cells expressing endogenous wild type or mutant EIF1AX. Ribosome occupancy of the global translational apparatus was sensitive to suppression of wild type but not mutant EIF1AX. Together, these studies suggest that cells expressing mutant EIF1AX may exhibit aberrant translational regulation, which may provide clonal selective advantage in the subset of uveal melanoma that harbors this mutation.

Drug Repositioning through Systematic Mining of Gene Coexpression Networks in Cancer.

Gene coexpression network analysis is a powerful "data-driven" approach essential for understanding cancer biology and mechanisms of tumor development. Yet, despite the completion of thousands of studies on cancer gene expression, there have been few attempts to normalize and integrate co-expression data from scattered sources in a concise "meta-analysis" framework. We generated such a resource by exploring gene coexpression networks in 82 microarray datasets from 9 major human cancer types. The analysis was conducted using an elaborate weighted gene coexpression network (WGCNA) methodology and identified over 3,000 robust gene coexpression modules. The modules covered a range of known tumor features, such as proliferation, extracellular matrix remodeling, hypoxia, inflammation, angiogenesis, tumor differentiation programs, specific signaling pathways, genomic alterations, and biomarkers of individual tumor subtypes. To prioritize genes with respect to those tumor features, we ranked genes within each module by connectivity, leading to identification of module-specific functionally prominent hub genes. To showcase the utility of this network information, we positioned known cancer drug targets within the coexpression networks and predicted that Anakinra, an anti-rheumatoid therapeutic agent, may be promising for development in colorectal cancer. We offer a comprehensive, normalized and well documented collection of >3000 gene coexpression modules in a variety of cancers as a rich data resource to facilitate further progress in cancer research.